Bicocca-FT-01-05 
hep-ph/01050T5 
May 2001 



Power Corrections in Flavour- Singlet 
Deep Inelastic Scattering 



G.E. SmyeQ 

Dipartimento di Fisica, Universita di Milano-Bicocca, 
and INFN Sezione di Milano, Italy. 

Abstract 

We investigate the 1/Q 2 power-suppressed corrections to structure functions in 
the flavour-singlet channel of deep inelastic lepton scattering arising from renor- 
malon insertions into an initial-state gluon, as obtained using the dispersive 
approach. The pinch-technique is used as a convenient tool in the separation 
of contributions. 
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1 Introduction 



We have learned much about the strong interaction of particle physics from deep inelastic 
scattering (DIS) experiments over many years. This process is continuing, with present and 
future experiments generating more data over larger regions of phase space. 

Alongside this ongoing experimental work, theoretical developments continue to be 
made. QCD, the established theory of strong interactions, continues to pose challenges: 
perturbation theory has been relatively successful at high energies, but the series expansion 
even here is at best asymptotic. In the non-perturbative regime lattice techniques are 
making advances. Yet a clean distinction between 'perturbative' and 'non-perturbative' 
cannot be made: all QCD observables involve some interplay between them. 

Although the structure functions (for example) cannot be calculated using perturbative 
QCD, the general shape of their asymptotic Q 2 behaviour is well known: the observed 
Bjorken scaling at high Q 2 is violated by additional smaller terms. The dominant scaling 
violation is a logarithmic Q 2 dependence, originating in the scale dependence of the parton 
density functions of the incident hadron. This can be calculated using perturbative QCD 
and used to measure the strong coupling a s . 

In addition to the logarithmic scaling violations, there are known to be contributions 
behaving as inverse powers of the hard scale, i.e. as 1/Q n - These include both correc- 
tions due to the non-zero hadron mass M, which are suppressed by a factor M 2 /Q 2 , and 
non-perturbative power-suppressed terms arising from higher-twist operators in the oper- 
ator product expansion. Such contributions are not included in fixed-order perturbative 
calculations, yet they are known to be important over the wide Q 2 range of available data. 

Over a number of years the higher-twist power-suppressed terms of a wide variety of 
observables have been estimated using two related approaches, the 'renormalon' (see for 
a review) and 'dispersive' [Q] models, although it is the dispersive approach that is used here. 
If we consider graphs with an arbitrary number of loop insertions in a gluon propagator, 
we assume that we can reconstruct a well-defined effective strong coupling at the scale of 
the gluon virtuality. The difference between this 'true' coupling and that reconstructed 
using fixed order perturbation theory gives rise to non-perturbative corrections, which are 
typically power-behaved. For a given observable, the shape of the leading correction can be 
found; an estimation of its magnitude requires the additional assumption of universality. 
Thus, starting from perturbative QCD, we aim to investigate the transition to the non- 
perturbative region. 

These approaches have been applied to various QCD observables, with the assumption 
of universality approximately holding ||. In flavour non-singlet DIS there are results for 
structure functions fragmentation functions and event shape variables |7[], while 

studies have been made of power corrections to structure functions ||, |§ and fragmentation 
functions [ID in the flavour singlet channel. 

The flavour singlet contribution to DIS is that involving the interaction of the gluons 
within the hadron. (We do not study the fermionic singlet, arising from the total sum of 
quark distributions, since these behave just like the non-singlet contribution j|, |5].) Power 
corrections are calculated using the renormalon or dispersive models by making insertions 
into the gluon propagators as shown in figure [I]. Here the lower part of the diagram 
represents the finding of a virtual gluon within the proton, and the upper part represents 
photon-gluon fusion. 

In the normal perturbative treatment of DIS, the asymptotic freedom of QCD enables 
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Figure 1: Power corrections in flavour singlet DIS 



us to treat the initial-state partons as free particles confined within the nucleon; so in a 
singlet calculation we would start with a free gluon and convolute the perturbative result 
with the gluon distribution function g(x). We cannot however do this in a calculation of 
power corrections, since the models we use consider modifications to the gluon propagator 
(loop insertions in the renormalon model, or, equivalently, a 'mass' in the dispersive ap- 
proach). We therefore consider the initial-state gluon to be generated by some perturbative 
mechanism, the simplest of which is by radiation from an on-shell parton. 

We immediately encounter a very serious problem: if the gluon is not on shell, both 
the upper and lower parts of this diagram are gauge-dependent. Gauge independence is 
achieved only when we include all the elements of some subset of diagrams of a given order 
in a s , and this includes diagrams that cannot be separated into the two halves of figure |l| 
diagrams which do not have analogous gluon propagators in which to make insertions or 
from which to take a scale for the running coupling. 

This is precisely the same problem encountered in the two-loop calculations of power 
corrections to event shapes in e + e~ annihilation, as discussed in fQ |, except that the 



gluon is now in a different channel. We might thus expect a similar solution: namely, 
the use of the pinch technique [0 to generate the running coupling at the scale of the 
gluon virtuality [13] in squared diagrams with two exchanged gluons, and the remainder of 
diagrams contributing to the power correction via some as-yet unknown mechanism. 

The first studies [§, [| [10|] of these non-singlet quantities considered the incoming virtual 
gluon to be radiated from a quark line, as in figure |2|. This is the simplest case, since the 
lower half of the diagram is gauge independent, as is the upper half to leading order, 0(a s ). 

There is however another problem to be overcome, concerning the interpretation of the 
lower half of the diagram. We may try to recover the singlet contribution to the power 
corrections by deconvoluting the full result with the quark-to-gluon splitting function, (this 
is the approach taken in ||), or we may leave the result as it is and interpret it as a genuine 
second-order contribution. These two interpretations give very different predictions for the 



2 



Figure 2: Radiation of gluon from quark line 



magnitudes of the power-suppressed corrections. I argue below (in section |^) that the latter 
approach should be adopted. 

In section |2| we discuss the flavour singlet contribution to deep inelastic scattering. The 
dispersive approach to power corrections is then briefly reviewed in section |3|. Section f| 
examines in detail the virtual gluon production represented by the bottom half of figure 
|l|, and the application of the pinch technique to this, while in section [| we calculate the 
contribution to the structure functions from photon-gluon fusion. Section || investigates 
the distribution of virtual gluons in the proton, paying particular attention to the 1/Q 2 
power corrections. Finally a summary is given in section [7[ 



2 Flavour singlet DIS 

Consider the deep inelastic scattering of a lepton with 4-momentum / from a hadron with 
4-momentum P, as shown in figure [| If the momentum transfer is q, the usual kinematic 
variables are Q 2 = —q 2 , the Bjorken variable x = Q 2 /2P ■ q, and y = P ■ qj P ■ I ~ Q 2 /xs, 
where s is the square of the energy in the cm. frame. 

Then the differential cross section is 

^ = ^{[l + (l-ym(x) + 2(l- y )F L (x)} , (2.1) 

where Fj-(x) = 2F 1 (x) and Fi(x) = F 2 (x)/x — 2Fi(x) are the transverse and longitudinal 
structure functions, which also have a weak Q 2 dependence which we do not show explicitly. 
(For simplicity we are neglecting any contribution from weak interactions, i.e. Z° or W ± 
exchange.) 

We can consider the photon to interact with an asymptotically free parton, moving 
collinearly with the hadron with momentum = xP^/l; (x < £ < 1). Then in the parton 
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Figure 3: Flavour singlet contribution to deep inelastic scattering 



model, to order a°, we have 

F T (x) = Y,e 2 q [q(x) + q(x)] F L (x) = , (2.2) 

q 

where q(x) and q(x) are the quark and antiquark distributions in the target hadron. Thus 
at this level there is no contribution arising from the gluon distribution in the hadron. 

The 0(a s ) contributions are most easily given as the distribution in the final-state 
variable r\ = P ■ r/P ■ q, (0 < r\ < 1): 

^ F i(x) = ^Y,4 £j{CFCit(t,v)[q{x/t) + q{x/$]+W , (2.3) 

where g(x) is the gluon distribution in the target hadron, Cp = 4/3, Tr = 1/2 and the 
coefficient functions Ci t j(£,r)) fl4j] are 

ClM>v) = ^r/ (2.5) 
Qr^v) = r + a-O 2 ] " 2 ^ 1 /^ 2 (2-6) 

c L>g {t;,v) = 8^(1-0- (2-7) 

An integration over the entire range of r\ from to 1 requires the implementation of a 
factorisation scheme to regulate the collinear divergences. In our calculations this will be 
effected by the introduction of a small gluon 'mass'. In addition the coefficient functions 
Ci >q (£, rf) acquire contributions at r\ — 1 from virtual gluon emission. 

We are interested in the singlet contributions Cj jS (£, rj), corresponding to photon-gluon 
fusion; the l/Q 2 power corrections from the non-singlet components Cj ig (£, r/) were success- 
fully analysed in 0, [j| . 
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3 The dispersive approach to power corrections 



We assume that the QCD running coupling a s (k 2 ) can be defined for all positive k 2 , 
and that apart from a branch cut along the negative real axis there are no singularities in 
the complex plane. It follows that we may write the formal dispersion relation: 

a s (k 2 ) = - f-^-p^), (3.1) 

JO /LT + k z 

where the 'spectral function' p s represents the discontinuity across the cut: 

Ps(/i 2 ) = ^-.ia 8 (ji 2 ^) - a s (/iV^)) = -^Disc a s (-p 2 ) . (3.2) 
Ztii L J Zm 

To lowest order in perturbation theory we have p s (p 2 ) — — ^\ a s(— P 2 )| 2 - 

Non-perturbative effects at long distances are expected to give rise to a non-perturbative 
modification to the perturbatively-calculated strong coupling at low scales, 5a s (p 2 ) = 
a s (p 2 ) — af T (/i 2 ), where af T (/i 2 ) is the perturbatively-calculated running coupling |fj]. 
Note that here a^ T (p 2 ) refers only to the contribution to the running coupling from a fixed 
(next-to-leading) order perturbative calculation, and so is itself well-behaved down to low 
scales, without any divergences: the Landau pole appears when we include an arbitrary 
number of loop insertions in the propagator. Hence Sa s (p 2 ) is assumed to be well-defined 
for all positive p? . 

We now consider the calculation of some observable F in an improved approximation 
which takes into account fixed-order contributions plus those higher-order terms that lead 
to the running of a s . As is well documented, for processes involving a single gluon, it 
is required to calculate the relevant contributions as though the gluon had a small mass 
fi 2 = eQ 2 . Calculations of power corrections in e + e~ annihilation, in non-singlet DIS, and 
in the Drell-Yan process, use a single such gluon. In these cases the 1/Q n corrections are 
found to be proportional to 



C F r°° d/j 2 
2tt Jo p 2 



^» = 75T /. J T^ n 5a s ^ 2 ). (3.3) 



Numerical values for these parameters must be obtained from data: fits of the 1/Q 2 cor- 
rections to DIS structure functions suggest that A2 — 0.2 GeV 2 [|J. 

However, DIS in the singlet channel involves two such gluons. Both gluons have an 
associated dispersive variable, so we obtain a characteristic function jF(ei,62), where e» = 
p 2 /Q 2 , which is simply the observable F calculated as though the gluons had masses \i\ 
and fi2, and without any factors of a s . 

Since both gluons are constrained to have the same 4- momentum s, we can simplify 
this to require only one dispersive variable. By defining p = —s 2 /Q 2 , we see that the 
dependence of the characteristic function T on e% and €2 is given by 

J (p + d)(p + e2) 

where the integration limits and the function / depend on the particular calculation. This 
may be re-expressed in the form 

F{ei,E 2 ) = , (3.5) 

ei - e 2 
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where 

So it is sufficient to perform the calculation with one 'mass' set equal to zero, the other 
giving us the form of the characteristic function. 
The alternative expression 

fi ic)= l -f dpm - 1 -f d -i£M ,3.7) 

e J e J p + e 

shows how this relates to the slightly different definition of T found in || IH|. In these refer- 
ences it is shown that a characteristic function of the form ( |3.5|) leads to a non-perturbative 
correction to the observable F given by 



lo n 

where, using Cauchy's theorem 



5F = | o °° % (2a s {fi 2 )8a s (fi 2 ) - [5a s (fi 2 )] 2 ) G^/Q 2 ) , (3.8) 



G(e) = -^DisoF(-e) = -JLj^e*) - f{ee^ )} = /(e) . (3.9) 

(Another method of arriving at (|3.8|) in this particular case is simply to insert the running 
coupling a s (—s 2 ) into the matrix element before integrating.) 

Since 8a s (n 2 ) is small in the perturbative regime, and vanishes as /i 2 — > oo, the correc- 
tion (|3.8|) depends on the behaviour of T at small e. So we perform an expansion of this 
function about e = 0. 

Any divergent term in T[e) is subtracted off. These are the terms causing the running 
of parton distributions, which gives rise to logarithmic scaling violations: the gluon mass 
here behaves as a regulator. 

The remaining terms, which vanish as e — > 0, give power corrections. We will find the 
dominant terms to be: 



J"~aieloge 5F = aj^r , (3.10) 

Q 



and 

f ~ |a 2 elog 2 e = 5F = a 2 ^| log ^| , (3.11) 



where the non-perturbative parameters D\ and D2 are defined by: 

dfj, 



D 1 = r%» 2 (2a s (n 2 )5a s (fi 2 ) - i5a s (n 2 )] 2 ) , (3.12) 

J fi 



\ogD 2 = T^^logfi 2 (2a s (fi 2 )5a s (fi 2 )-[5a s (fi 2 )] 2 ) . (3.13) 

D\ JO j2 z v ' 

While we expect the form of a s (// 2 ), and hence and D 2 , to be universal, we have as 
yet no numerical values for them. It will be necessary therefore to extract values for Di 
and D 2 , either from experimental results or from some model of the form of a s (fi 2 ). 
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4 Production of the gluon 



4.1 Gluon kinematics and polarisations 

Let us consider the radiation of a virtual gluon with 4-momentum s from a massless 
parton with 4-momentum p, as shown for the case of a quark in figure |2|. The parton 
model assumption is p = xP/£, where P is the 4-momentum of the incoming hadron. Since 
the gluon is virtual, it need not be collinear with the proton, nor need it be transversely 
polarised. 

We describe the kinematics by three additional variables, u and £' defined by 

e q 2 

Ta=-7^=P> (4-1) 



l+uj 2s -q £' 2 ~~ Q 2 

and 9 the azimuthal angle in the Breit frame. Then the relevant 4-momenta are, in the 
Breit frame, 

q = ±Q(0,0,0,2) (4.2) 

p = ±Q(l/£,0,0,-l/£) (4.3) 

s = ^Q(s , s±cos9, s±sm9, s 3 ) , (4.4) 

where 

so _ i ± (i - m> m 

4 = ^ziim^m (4 . 6) 

■H = (4.7) 

The gluon is therefore produced with three degrees of freedom: £' determines the longi- 
tudinal momentum fraction x/£' of the gluon, u is proportional to the gluon's virtuality, 
and 9 gives its azimuthal angle. If we are concerned only with the produced gluon, and 
not with the parton that emitted it, then £ is also free. Note that in the limit u — > 0, 
the gluon 4-momentum s becomes xP/^ ; , the usual parton model result, independent of £ 
and 9. This limit gives the standard photon-gluon fusion contribution to the observable 
under consideration. However in section |3] we saw that it is terms of higher order in p that 
generate power corrections: the transverse momentum of the gluon is crucially important. 
Consequently we might like the gluon to have a 4-dimensional 4-momentum distribution 
within the proton: this will be discussed in section ||. 

With these definitions, the integration measure for the splitting becomes 

d 3 p' 1 f 1 de' f 2 * d9 * du 2 £ 

= WfhYh Y,L - (1 " } e • (4 - 8) 

Let us also introduce polarisation vectors ef, (i = 1,2,3), for the radiated gluon. We 
must choose the ef to satisfy s ■ = and £j • €j = — 5y. Then we have the identity 

£« = -<r + £ 4- • (4-9) 
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For those contributions to figure [l] where both the top and bottom halves are gauge- 
independent (such as the lowest-order contribution to figure 0), we can use ( |Q| ) to replace 
the gluon propagators with sums over polarisation vectors, thereby detaching the top and 
bottom halves of the diagram. Where it is not gauge independent, this gives the contribu- 
tion in Landau gauge. 

For the purpose of later discussion, let us introduce two such sets of polarisation vectors. 
Firstly, let us define the 'natural' polarisation vectors ^ by imposing g-e 2 = q-e 3 = p-e 3 = 0. 
Then the identity ( |L9|) gives us 



iQ(l - Uj)y/U 



£ m = [(1 + uj)C - 2u£\ s» + [2 - (1 + u^/e\u g»-(l- co)^p^ 



Q(l-u)yJu(l-Z/?)(l-u£/?) 

and €3 is orthogonal to the vector space spanned by s, q and p. 

The vectors e 2 and 63 are both orthogonal to q, and so are the transverse polarisation 
vectors in frames in which s and q are (anti-)parallel, and thus in all physical frames. Also, 
since s is spacelike, the longitudinal polarisation vector e% is timelike. (The factor % in the 
expression for e\ then causes the identity Q4.9|) to hold as written; we could alternatively 
have introduced a factor —1 in the sum over polarisation vectors.) 

These are the natural polarisation vectors because e\ is longitudinal and e 2 and €3 are 
transverse: they are thus the polarisation vectors that should be used when considering 
parton splitting, and the splitting functions arise naturally from their use. 

Secondly, let us define the 'diagonal' basis of polarisation vectors e, by imposing p ■ e 2 = 
p ■ e 3 = q ■ e 3 = 0. It will be seen below that this is the basis that diagonalises the matrix 
describing gluon production, and is formed by mixing and e 2 to give the new basis 
vectors 
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f 1 ' 

f 2 



iQy/uj 

^'-(l+^ei^-e^-K/eo <t 



Q^u{l-t/t>){l-Lot/t>) 

~4 = 4- (4.11) 

The disadvantage of this basis is that <L\ and e 2 are neither transverse nor longitudinal, and 
therefore are not suitable for a discussion of parton splitting or of splitting functions. The 
use of this basis will be denoted by a tilde on appropriate symbols. 



4.2 Radiation from quark line 



The matrix element for the process shown in figure integrated using the measure 
jf), gives the following parton- level contribution to some observable F: 



2ir 



'■j 



1 d£' /** d9 ft' du 



5 £' Jo 2tc Jo uj 



2 % 3 



2{p • €i){p ■ €j) 



(4.12) 

where Fy is the contribution due to photon-gluon fusion with the appropriate gluon polar- 
isations. 
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Using the natural basis ( [4.10 ), we find 
where 



(4.13) 



( -2(l-z)(l-^z) (2-z-wz)y/(l-z)(l-wz) 

(1— £j)z i(l—ui)z 



M?Az,u) 



V 



(2-Z-Uz)y/ (1-Z)(1-WZ) 

i{l—u>)z 





{2-z-mz) 2 

2(l-£j)2 









(1-LJ)Z 

2 



(4.14) 



Several comments may be made about this matrix: 

(i) The elements on the diagonal correspond to the exchange of a gluon with definite 
polarisation; the two off-diagonal elements give rise to an interference term between 
the longitudinal and one of the transverse polarisations. 

(ii) The usual q — > g splitting function can be recovered from the total transverse piece 



Mb + Mas in the limit 



0: 



lim (Mf 2 + Mf 3 ) 



l + (l-z) 



(4.15) 

The standard factorised expression for the leading (logarithmic) divergence arises 
from the fact that as u — > 0, s M — > xP^/£' and ~ s^/y/s^, so F n , F 12 and F 2 \ 
all become subleading and (for an unpolarised observable) F 2 2 and F 33 become equal 
and functions of £' only. Thus the leading logarithmic divergence in F is 

^ leadms (o - ^ f f f -p^,(e/o^(0 (4.i6) 

Z7r j§ 4 "'o u; 

where C 9 represents the contribution from a transversely-polarised initial-state gluon. 

This however applies only to the leading divergence: other pieces, in particular those 
non-divergent pieces generating power corrections, do not factorise in this way. In 
general, the functions Fij have dependence on £ and 6 as well as uj and i.e. for the 
purposes of higher twist contributions the gluon 'remembers' the momentum fraction 
of the particle that emitted it. 

(iii) The matrix Mfj has zero determinant, and hence has a zero eigenvalue. The other 
two eigenvalues are both equal to (1 — u)z/2. It is diagonalised if we use the diagonal 
basis ( |4.11| ), in which case the decomposition ( 4.12|) simplifies considerably to: 

pq = OsCf r 1 d£_ r 2 * d0_ re du (1 - m )Z 
2tt 7{ f Jo 2vr Jo u 2f 



|f 2 2 + -F33I 



(4.17) 



In other words, only two gluon polarisations, e 2 and e 3 , are permitted, but they are 
not the transverse polarisations assumed in the parton model. Since this basis fails to 
make apparent the form of the leading divergence and does not naturally give rise to 
the splitting function, it is clear that any attempt to remove the initial quark line from 
the power corrections calculated in |5| and [ID] cannot be simply a matter, as with the 
logarithmic divergences, of deconvoluting with the quark to gluon splitting function. 
Also, since all the elements of ( |4. 14[ ) contribute to the power correction we would 
need to be able to include the unphysical polarisations in the gluon distribution. 
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4.3 Radiation from gluon line 



The radiation of the virtual gluon from a gluon line is shown in figure |j. This diagram 
is gauge-dependent, and while its logarithmically divergent piece is gauge- independent and 
given by the g — > g splitting function, the remainder of the diagram, and in particular the 
pieces generating power corrections, are not. 





Figure 4: Radiation of gluon from gluon line 



The matrix element for the process shown in figure [|, integrated using ( f4.8|) and using 
the natural basis of polarisation vectors ( ^.1U| ), gives the following contribution to the 
observable F: 



F = — — 2^ / TT / oZ / — M m/i > i , u, 

4 Jo z?r Jo w 



2tt 



'■j 



where 



2z 



(2-z)(2-2z+z 2 ) 
2izWl—z 



(2-z)(2-2z+z 2 ) 2-4z+4z 2 -2z A +z 4 
2iz^/l—z z) 














z(2-2z+z 2 ) 

1-2 



0(u>) 



(4.18) 



(4.19) 



To compare with the quark case, several comments should be made: 



(i) The terms of 0(u) are gauge-dependent, and therefore of not much use until we can 
isolate pieces of other diagrams to give a gauge-invariant total. The pinch technique 
(see next section) succeeds in giving a gauge- invariant result, but at the expense of 
changing this leading order matrix. 

(ii) Again there are on-diagonal elements corresponding to exchange of a gluon of definite 
polarisation, and off-diagonal elements giving interferences. In fact all six off-diagonal 
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elements are non-zero at 0(u) level, unlike the quark case. In addition the use of the 
basis ef does not diagonalise the matrix, even to leading order. 



(iii) The transverse elements M| 2 and M| 3 yield polarised splitting functions, and their 
sum gives 

lim (M| 2 + M| 3 ) = ^'^p 2 = P 9 -> 9 (z) • (4.20) 
Thus the leading divergence is the well-known result 

^ leadins (0 = ^ f f f - W£/03,(0 ■ (4-21) 

Z7T J£ £' Jo UJ 

(iv) The matrix Mfj diverges at z = 1, i.e. when £' = £. This is due to the gluon with 
momentum p' becoming soft, and the divergence is cancelled by virtual corrections. 



4.4 Use of the pinch technique 



When we considered the diagram in figure |], we found that the terms giving rise to 
power corrections (i.e. the terms in the integrand of ( J4.18| ) that do not diverge as u — » 0) are 
gauge-dependent. We must consequently include pieces of other diagrams to restore gauge 
invariance. We could add in all 0(a 2 ) diagrams that contribute to the process 7*g — > qqg, 
as in 'Milan factor' calculations for event shape variables (see [|11|), but that would force 
us to include diagrams unlike that of figure |l|: it is not clear what relation (if any) these 
diagrams have to the scale chosen for the running coupling. 

First we notice that figure [| is related by crossing to the e + e~ annihilation diagram 
shown in figure [S]. So, just as in e + e~ annihilation in we can apply the pinch technique 
to the diagram to define a gauge-invariant contribution. But while the pinch technique 
has previously been applied to internal gluon loops JT3] and cut loops containing outgoing 
partons ||11|| , it is now being applied to a cut loop where one of the particles is incoming 
and the other outgoing. This difference is irrelevant since the algebra is identical. 




The matrix element for figure on applying the pinch technique, is given in equation 
(A. 9) of [11]: in our notation it may be written 

|M| 2 = g 2 C A \ {8(s 2 g a p - s a sp) + 2(p a + p' a ){p p +p' /3 )}j2 , (4.22) 

S i,j 
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where F^ is the contribution from the shaded blob. Using the crossing relation and inte- 
grating using ( |4.8| ), we obtain the following form for the observable F: 

2(p • 6i)(p ■ Sj) 



2tt 



'•J 



1 dg f 2 - d9 ri' dw. N £ 

1 — UJ) 



£ £' Jo 2tt Jo uj 



2b ij 



Fait,?,", 



(4.23) 

Note that this differs from the incoming quark case only by having a different constant 
multiplying 5ij. In fact the only linearly independent terms allowed inside the square 
bracket, consistent with gauge invariance and having the correct dimensions, are the two 
shown. This was also seen in the e + e~ annihilation calculation of [[IT . 
Using the natural basis ( f4.10|) , we find 



2n 



1 d? f 2 - d9 /■€' duo 



£' Jo 2ir Jo uj 



where 



( -(2-3z+ljz)(2+z-3ljz) (2-z-uz)yJ (l-z)(l-uz) 
2(1— w)z i(l—u>)z 



MlAz,u) 








(2-z-uz)y/ (l-z){l-wz) (2-z-ujz) 2 _|_ 3(l-u))z 
i(l-u)z 2(l-u)z ' 2 

2(1 -uj)z j 



(4.24) 



(4.25) 



As previously, there are several comments to be made: 

(i) Although this matrix is gauge invariant, it does not reduce to ( 4.19Q in the limit 
uo — ► 0. One consequence of this is that the g — > g splitting function is no longer 
trivially recovered from the total transverse piece M| 2 + M| 3 in the limit uj — » 0: 



lim 



22 



+ MI 



2(1 - z + 2z 2 



33; 



^0 ■ - z 
where the usual splitting function is 

2(1 - z + z 2 ) 2 



2z 3 



1-z 



g-*g 



z(l-z) 



(4.26) 



(4.27) 



The additional piece arises from the fact that the pinch part of the diagram was 
removed (or alternatively that the pinch parts of other diagrams were added) in order 
to secure gauge invariance. This can be understood as follows: the factor 1/(1 — z) in 
the splitting function can only arise from the collinear limit of the gauge- dependent 
piece (n ■ p)/(n ■ p') in the sum over polarisations of the p 1 gluon in figure |j. In 
order to achieve gauge invariance we must add or remove terms so as to cancel all 
n-dependence, not just in the collinear limit but identically for all p and p', and so 
terms with the factor 1/(1 — z) must necessarily disappear. 

This happens even though the remaining diagrams themselves have no logarithmic 
divergence associated with the gluon splitting (although of course they may have 
divergences arising elsewhere): applying the pinch technique separates the diagrams 
into a pinch part and a remainder, which have equal and opposite logarithmically 
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divergent pieces. Thus this is not a sensible way to study splitting functions. Also, 
since the only two linearly-independent gauge-invariant terms with the correct di- 
mension are those in ( |4.23|) , and no linear combination of these can reproduce the 
splitting function, we see that we cannot have a fully gauge invariant expression for 
the splitting g — > g+g* (with some suitable definition for what is meant by g*, e.g. use 
of pinch technique) while retaining the full splitting function. 

We also notice that the contribution ( |4.26| ) is no longer invariant under z — > 1 — z: 
this is because the two daughters of the splitting are no longer identical, one being a 
real on-shell gluon and the other a virtual gluon with a modified propagator. 

(ii) The difference between the production matrix employing the pinch technique and 
that using standard perturbation theory is 

/ 2z d2^k \ 

2«vl— 2 



-(2-z)z -z 3 q 

2i^f\=z 1-2 

V o o=^y 

\ 1-2 / ij 



+ 0(u) (4.28) 



(where 0{uj) terms are gauge-dependent). This is the contribution that gives the ad- 
ditions to the splitting function. It becomes zero as z — > 0, which is as the interacting 
virtual gluon becomes soft. 

(iii) As in the quark and conventional gluon cases, there are both on- and off-diagonal 
elements. Now, however, the matrix Mfj is diagonalised using the diagonal basis 
( [4. 11| ), just as in the case where the emission is from a quark line. (Indeed one can 
trivially see that any linear combination of the two terms in ( 4.23| ) is diagonalised 



using this basis.) One eigenvalue is 3(1 — u)z/2, and the other two are both 2(1 — to)z. 
Then the decomposition ( |4.23|) simplifies considerably to: 



2vr Jt f Jo 2n Jo u f 

5 Structure functions 

Consider now the part of the process involving the interaction of the virtual gluon 
within the proton with the virtual photon. This is the process shown in diagram |3|, which 
corresponds to the quantities above. 

Working in the Breit frame, the 4-momentum of the incoming photon and gluon are 
given by ( |4 . 2| ) and ( |4.4j ); in this part of the process the variables uo and 9 are considered 
fixed. Let us introduce the variables T] = P ■ r/P ■ q, fj — P ■ k/ P ■ q , x the azimuthal angle 
between r and s, and x the azimuthal angle between k and s. The kinematics are then 
given by: 

lQ(z , z ± cos(x + 0), z ± sin(x + 9), z 3 ) (5.1) 



r 



2 



k = lQ(z ,z±cos(x + S),z±sin(x + 9),z 3 ) . (5.2) 

The definitions of rj and fj along with the on-shell conditions for the outgoing particles 
require that 

z 2 z 2 
Zo = r] +^ Z3 = 7] ~^ (5 ' 3) 
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whence conservation of the Oth and 3rd components of 4-momentum give the conditions 

n+n = (5.5) 
o-mm-j/e), (5 . 6) 



7 2 =2 
Z ± | Z ± 

Arj Arj 



while conservation of transverse momentum requires that s±, z± and z± satisfy the triangle 
inequalities 



\z± — Sj_\ < Zj_ 
\z± — s±\ < z± . 



(5.7) 
(5.8) 



An additional variable, f3, is required to parametrise the permitted values of z± and z±. 
Let us choose to write: 



where 



z]_ = 477 
z\ = Arj 



arj + brj + 2 cos Pyabqr) 
ar] + bfj — 2 cos /3 '\J 1 abrjfj 



(5.9) 
(5.10) 

(5.11) 



e /2 (i-^/e 2 ) 2 e 2 (i-^/e 2 ) 2 

Given zj_ and the angles x an d X are determined up to a sign. We may then 
choose i] and (3 as the independent variables, with phase space 

< v < 1 - ^/r 2 



< < 7T . 

To integrate the matrix elements we apply the operator 

' J " k (2n) 4 6 4 (q + s-r-k) 



I 



(27r) 3 2r° (27r) 3 2fc° 
Integrating with respect to k and making substitutions for r gives 

1 r z ± dz ± dxdr] S(A- y/B - 2s ± z ± cosx) 



(5.12) 



(5.13) 



4(2vr) 



77 



y/B — 2s±z± cosx 



(5.14) 



where A and B do not depend on x- 

Next we integrate over x- There are two values satisfying the integration condition, and 
they differ by a sign. This gives 



drj dz 2 j 



77 s±z± \ smx\ 



4(2tt) 2 

Applying the parametrisation in terms of /3, we find that 

dz 2 t 



d(3 



= 81] sin j3 \J abr/rj 
s±z±\sinx\ = 4(1 - uji/i' 2 ) sin p^J abr/rj , 
14 



(5.15) 

(5.16) 
(5.17) 



and therefore the integral operator for the photon-gluon fusion is 



8tt 2 (1 - ui/i' 2 ) Jo 
We calculate the quantities u) as defined by 

a s T R 



drj / dp 



2tt 



^ dd „ _ w 

Z7T 



(5.18) 



(5.19) 



where the sum over g' represents the outgoing quark- ant iquark flavours. We first perform 
the 9 integral: the contributions to Fl have no ^-dependence, while the contributions to Ft 
have terms proportional to cos# and cos 2 9. Then the integration over (3 can be performed 
by writing t = tan(/3/2), and that over rj by carefully expanding in powers of u = p£ /2 as 
described in M. We obtain, for the longitudinal structure function, 



ii 



C 



cl 2 = cf 
cf 



cf 



32(1 - i') [(3 + C - 4f') + (1 - £') \ogu]uj + £>(cu 2 ) 
-32iJl - 0(1 - 3C' - f'logwV + 0(^ 2 ) 



8^(1 - + 8 (!/£' + 2f /£' - 10 - 8£ + 23£' + 9££' - 17£' 2 ) 



+(£/£' - 3 - 2£ + 8£' + 2££' - 6£ /Z ) logwju; + 0(cu 2 ) 

8£'(1 - £') + 8 [(!/£' + 2£/£' - 10 - 8£ + 23£' + 7££' - 15£' 2 



+(£/£' - 3 - 2£ + 8£' + 2££' - 6£' 2 ) log^cu + 0(^ 2 



Hence in the limit u — > 
Cf = Cf = 8^(1-0- 



we recover the well-known results Cf 



Cf 



Cf 



(5.20) 
(5.21) 

(5.22) 

(5.23) 
and 



For the transverse structure function we obtain 

-8£'(1 - £') + 8[(-2/3£' + 25/3 + 2£ - 15£' - 2££' + 8£' 2 
+2(1 -Cf\oguj\uj + 0{u 2 ) 



/-til 



(5.24) 



^12 W21 
0^ — 



^22 



-4i 



£/£' 



(£/3£' 2 - 2/3£' - ll£/3£' + 10/3 + 18£ - 20£' 



-14££' + 16£' 2 ) - 2(1 - - 3£) log u]u + 0(u? 



(5.25) 



-2(l-2£' + 2£ ,2 )(2 + log^)- 



(£ 2 /3£' 3 + 8£/3£ /2 - 8/3f 



+ 13£ 2 /3^ - 88£/3£' + 76/3 - 14£ 2 /£' + 60£ - 44£' + 16£ 2 - 48££' 
+30£ /2 ) + 2(£/£' 2 - 1/0 + i 2 lf - 7£/r + 6 - £ 2 /£' + 9£ - 7£' 

~/2 



^33 



+£ 2 - 6££' + 40) logcujcu + 0(u 2 
-2(l-2£' + 2£ /2 )(2 + log^) 



(5.26) 



1 - £/£' 



£ 2 /3£' 3 + 4£/£' 2 -4/£' 



/2 



+ ll£ 2 /3£' - 28£/£' + 24 - 18£ 2 /£' + 68£ - 52£' + 16£ 2 - 48££' + 34£ 



/2x 



+2(£/£" - + £ 2 /£" - 7f /£' + 6 - 3£ 2 /£' + 13£ - llf 
+3£ 2 - 10££' + 8£' 2 ) log u]u + 0(u 2 ) . 



(5.27) 



As u; — > the pieces C T 2 and Cf. 3 diverge logarithmically: this is the piece associated with 
the collinear splitting of the gluon into a quark-antiquark pair. 
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6 Power corrections 



We noted above that the interaction of the produced virtual gluon is given independently 
of the production mechanism by the function Fij(£, uo, 9), where i and j represent polar- 
isations and £, to and 9 parametrise the gluon's 4-momentum. In general, calculations 
of power corrections will depend on all four of these quantities. 

Thus a natural and intuitive way to combine this with a gluon distribution function 
might be to write 

n*) = £ f 1 d 4 t d 4 rf- f - ^/i'^iimA^'^) , (6.i) 

77 Jx f if C JO Ztt Jo uo 



hJ 



where gij(x /£',&, £/£') is the relevant gluon distribution function which by symmetry does 
not depend on 9. Of course this needs to be treated with some care, since both halves of 
the diagram in figure [l] are gauge-dependent and so therefore are all these distributions 
when the gluon is off-shell, i.e. w ^ 0. Nevertheless, we can still make progress: we present 
an intuitive argument that, although mathematically not totally rigorous, is helpful in 
understanding the underlying physics. 

For the diagonal elements we can view this probabilistically: the expected number 
of gluons of polarisation i in an element of parameter space duduodz at (u, uo,z), where 
u = x/£' is the longitudinal momentum fraction of the proton carried by the gluon, and 
z = 

gu(u,u,z)du . (6.2) 

00 z 

This is intended to be schematic onlyQ — clearly any gluon not collinear with the proton, 
i.e. with uo ^ 0, will experience the confining effect of the QCD potential, and so cannot 
exist for more than a short time (where 'short' in this context means ~ 1/Aqcd)- Such 
particles are not asymptotically free: they can only be resolved at high momentum scales, 
and it is precisely this resolution of gluons with non-zero transverse momentum that gives 
rise to the running of the parton distributions. Power corrections are also known to arise 
from the interactions of gluons with non-zero transverse momentum (see e.g. ||15||), and so 
are also generated by the resolution of short-lived virtual gluons within the proton. Thus 
we expect power corrections in singlet DIS to be intimately related to the running of parton 
distributions. 

Consider first the transverse polarisations: these are the physical polarisations. Speak- 
ing schematically (i.e. not being too precise about our choice of factorisation scheme), we 
may view the polarised gluon distribution function gi(u,M 2 ) as representing the gluons 
with transverse momentum less than M. So 

/ „ ,2n r 1 dz r (M2) duJ i 

gi{u,M ) = / — / — gu{u,u,z) , (6.3) 

Ju Z Jo UO 

for % = 2, 3, where uo(M 2 ) is the value of uo corresponding to a transverse momentum M 2 . 
Therefore, 

, /2 %M5 f^dzM 2 du 

M ^771 = / —T7Tj^gu{u,uo,z) , (6.4) 



DM 2 Ju z uo dM 2 

2 It could of course be more rigorously formulated in terms of well-defined unintegrated parton distribu- 
tions, but that is unnecessary for the purposes of the argument given here. 
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and, using equation (|4.6|) , we obtain 

M 2 duo 



1 — UOZ 



uo dM 2 \-2wz 
Further, we know from the DGLAP equation that 

a, r 1 dz r 



1 + O{uo) 



M2 d 9i (u,M 2 



dM 2 



2n Ju z 



C F q tot (u/z)P^ g (z) + C A g(u/z)P i 



(6.5) 



(6.6) 



so for the sake of illustration let us make the simplest consistent assignment for the differ- 
ential gluon distribution ga(u,uo,z), which is 



gu(u,uj,z) = — 

Z7T 



C F q tot (u)P^ g (z) + C A g(u)P^ Q (z)\ + O{uo) . 



(6.7) 



Here gtot(w) = SJ<?( M ) + q{ u )] is the total quark and antiquark content of the proton. 
So let us return to (|6.1| ) and define the coefficient function 

C g {Z')=F H (U'A9) (6.8) 

which is independent of £, 9 and polarisation i. Let us separate off the piece that generates 
the collinear divergence by writing 



de 



/ F H (Z,?,u,0)^- = C g (Oe(uo(M 2 ) -co) + 8F ii (£,£',u) 

JO 271 



(6.9) 



we then obtain the contribution to the observable F from incoming transversely-polarised 
gluons as: 



^'"trails (■£) 



E 

i=2 



d£, /-i d? /•«' duo d9 



x £ Je, £' Jo uo Jo 2tt 



gu(x/g,u,Z/g)F u (Z,Z',u, 



(6.10) 



1 dP 

-j r g{x/i\M 2 )C g {0 + 

x 5 



E 

i=2 



d£, f 1 d? f e duo 



i k e h uo 



gil (x/i\uo^/P)8F u ^i\uo) . 



(6.11) 



The first term on the right hand side of ( 6.1 lj ) is simply the standard contribution to 
photon-gluon fusion, evaluated with an on-shell initial gluon and convoluted with the gluon 
distribution g(x) = J2^=2 9i( x )- This contains logarithmic scaling violations given by the 
running of the gluon distribution, but contains no power corrections associated with the 
virtuality of this gluon. 

The second term gives us the required power corrections, since this is the term with an 
integral over the gluon's virtuality. Let us now substitute the differential gluon distribution 
with expression fl6.7Q, thus obtaining 



8F tiaT1B (x) = J2 I ~f f -7T f —9ii{.x/i',uo,i/i')5F li {i,i\uo) 
~ZJx £ J£ £ Jo v 



i=2 
3 ,1 



2ixj 



E 



i=2 



x 



S k i 

d£, [ l di' rt'dwf 

jj £ yio ~uo~ 



C F q tot {x/£,)P*MIO + C A g{x/£)P'{ £./£,') + O(uo 



(6.12) 



SFu. 
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Although the above argument has been restricted to the contributions from transverse 
polarisations, and even these have been treated only schematically, we are now in a position 
to be able to see what is happening physically We have gone round in a big circle: starting 
with figures |2] and |, we have detached the lower parts of the diagrams corresponding 
to gluon production, to give the piece corresponding to photon-gluon fusion. But then 
in order to convolute with the correct differential gluon distribution within the proton 
we had to make use of the DGLAP equation, which reintroduced those lower legs of the 
diagrams. Thus for the purposes of power corrections we find that the lower parts of those 
diagrams genuinely are important and cannot be removed. This also explains the presence 
of the longitudinal gluon polarisation: we are not dealing with a real incoming gluon and 
so we need not be restricted to physical polarisations. Equation ( |6.12| ) is intended to be 
illustrative only — it was not rigorously derived. Yet it is clear that while the leading 
perturbative contribution to singlet DIS, given by the first term in equation ( |6.11| ), is 
0(a s ), the 1/Q 2 power corrections as well as the logarithmic scaling violations are (9(a 2 ), 
the additional factor of a s arising from the DGLAP equation. (Contrast this with flavour 
non-singlet DIS, whose lowest order is O{a s ) but which has leading power corrections at 
0{a s ).) 

Thus the leading singlet power corrections are from the diagrams |2| and |] as they appear: 
we have two independent renormalon chains, and the magnitude of the power corrections 
is given by the quantities D\ and D 2 of equations ( |3.12| ) and ( |3.13|) . 



The characteristic function J-"(e) generating power corrections to structure functions is 
thus given by 

H*) = IE 4 I t J {^fltot(x/OC^e) + , (6.13) 




where 



C^^e) = f -^ r Y,^l 9 ^/^pf)C ij {^',pf) , (6-14) 



and the C 1 - 7 were defined in equation ( |5.19| ) and calculated in ( |5.20| ) to ( |5.27| ). The sum 



over q' represents outgoing quark- ant iquark flavours in the photon-gluon fusion. We retain 
terms in J^(e) up to 0(e) that are non-analytic as e — ► 0. In order to define a gauge- 
invariant quantity with a natural scale for the running coupling, we use the matrix (|4.25| ) 



evaluated using the pinch technique — the remaining diagrams are also expected to give 
power corrections, but as in (Tl| they are expected to have a different structure. 
Using ([CTBp and (fO|) we find * 

^C q T = -|(2-63e + 63^-2^ + 121oge-27^1oge-27^1oge + 12^ 3 log01oge 
+| (4 + 3£ - 3£ 2 - 4£ 3 + 6£ log f + 6£ 2 log £) (log £ - 1 + § log e) log e 
+ |(2 + 25£ 2 - 25£ 3 - 2£ 5 + 15£ 2 log£ + 15£ 3 logfjeloge 

-2(5£ 2 - 5£ 3 + 2£ 2 log£ + 2£ 3 logOQog£ - 1 + |loge)eloge , (6.15) 
iCl = -§(l-3£ + 2£ 3 -3£ 2 log01oge 

-2§s(17 + 75£ 2 - 125£ 3 + 33£ 5 + 301og£ + 75£ 3 log£ - 45£ 5 logOeloge 
+^(2 - 15£ 2 + 10£ 3 + 3£ 5 - 15£ 3 logO(log£ - 1 + iloge)eloge , (6.16) 

which is the result given in ||, and 

= -|(l-99£ + 99£ 2 -£ 3 + 61og£-27£loge-54e 2 loge + 33£ 3 log01oge 
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+ |(2 + 6£ + 3£ 2 - lie 3 + 3£log£ + 12£ 2 logO(k)g£ - 1 + A loge) loge 
+ f(l + 40£ 2 - 35£ 3 - 6£ 5 + 15£ 2 log£ + 30£ 3 logOdoge 
-8(5£ 2 - 5£ 3 + 2£ 2 log£ + 2e 3 logO(loge - 1 + |loge)eloge , (6.17) 
£C 9 L = -|(1 - 3£-9£ 2 + 11^ -12£ 2 k)gO loge 

+750£ 3 log f - 270£ 5 log £)e log e 
+Aj(l - 15£ 2 + 5£ 3 + 9£ 5 - 30eMogO(log£ - 1 + iloge)eloge . (6.18) 

Thus we obtain the 1/Q 2 power corrections to Ft(x) and Fl(x): those arising from the 
quark distribution are 



where 

5C q T (0 



Di 

Q 2 



Di 

Q 2 



2(2 + 25£ 2 - 25£ 3 - ie + 15£ 2 log£ + 15e 3 logO 



(6.19) 



5£ 



-2(5f - 5£ 2 + 2£log£ + 2£ 2 logO log 



F>ij 

eQ 2 



(6.20) 



8(17 + 75£ 2 - 125£ 3 + 33£ 5 + 30 log£ + 75£ 3 log£ - 45^ 5 logf) 

225e 

lo § -7^ 



and those from the gluon distribution are 

TrCa 



x 



eQ 2 



g{x/05C!{(i) 



(6.21) 



where 



(27T) 2 \^ * I J, L 
4(1 + 40£ 2 - 35£ 3 - 6£ 5 + 15£ 2 log f + 30£ 3 log 



(6.22) 



5£ 



-8(5f - 5£ 2 + 2£ log £ + 2£ 2 log log 



D 2 i 



eQ 2 



8_ m + 1275£ 2 - 1475£ 3 + 183£ 5 + 30 log£ + 
Q 2 L 225£ V450£ 2 log f + 750£ 3 log £ - 270£ 5 log £ 
, 16(1 - 15£ 2 + 5£ 3 + 9£ 5 - 30e 3 logQ . D 2 £ 

+ — log -Tyl 



15^ 



eQ 2 



(6.23) 



(6.24) 



7 Results and conclusions 

Figure |] shows the magnitudes K of the 1/Q 2 power corrections to structure functions, 
given by 

5F«%{x) _ Di 



,(0), - n2 K T/L( X ) ' 



4 0) (x) Q 1 



(7.1) 



19 



where F T °\x) is the Born-level result for the transverse structure function given in equation 
fl2.2|). These were calculated at Q 2 = 500 GeV 2 , using the corresponding MRST (central 
gluon) parton distributions ||16|| , and assuming four flavours of outgoing quark- ant iquark 
pairs. The unknown value of D 2 /e was set to 0.06 GeV 2 , i.e. approximately A 2 , although 
the qualitative behaviour of these power corrections does not change provided we keep 
D 2 < Q 2 . 



-l 



-3 



V K T (x) 




K L (x) 



0.05 



0.1 



0.15 x 0.2 



0.25 



0.3 



0.35 



Figure 6: Graph showing Kt{x) and Kl{x). The dotted lines represent the contributions 
from the quark distribution; the dashed lines represent the gluon contributions. The totals 
are represented by the solid lines. 



The dominant contribution to both of these comes from the gluon distribution, as may 
be expected at these relatively low values of x. All the contributions tend to zero at large 
x, and become large at small x. The contributions to K T are positive and those to K L 
negative, but they have similar magnitudes. 

The corresponding quantity related to F 2 /x, which is K 2 /x = Kt + K^, is shown in 
figure [7|. The positive and negative contributions partially cancel, giving a power correction 
smaller by a factor 3 or 4 than that for Ft or Fl- 

The power corrections are multiplied by the unknown factor D\, defined in ( 3.12j) . If 
Di is positive, we might expect Ft(x) to show a positive 1/Q 2 power correction and Fl(x) 
a negative one. These results are all qualitatively the same as those in |S]], which took into 
account only the quark contribution, but the inclusion of the gluon contribution increases 
considerably the size of the correction, in the case of Ft by more than an order of magnitude. 



8 Discussion 



The calculation of power corrections using the renormalon model or the dispersive 
approach in flavour singlet DIS is a non-trivial problem, due to the fact that figure [I] as 
it stands is gauge dependent, both in the upper and lower halves of the diagram. These 
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Figure 7: Graph showing Kz{x)/x. The dotted line represents the contribution from the 
quark distribution; the dashed line represents the gluon contribution. The total is repre- 
sented by the solid line. 



halves may both be considered to be cut insertions into a single gluon propagator, albeit 
in another channel, and so the pinch technique can be used to define a natural gauge- 
invariant quantity. The diagram [I] evaluated using the pinch technique gives only one 
gauge-independent part of the full 0(a 2 ) amplitude, and the expectation is that the other 
parts will also give power corrections, although not taking the form of two renormalon 
chains with equal 4-momenta. 

Power corrections thus arise when a virtual gluon is radiated from a quark or gluon, and 
this virtual gluon then interacts with the photon. We found that, as well as the physical 
transverse polarisations, we also have contributions to power corrections from longitudinal 
polarisations and interference terms. In addition, since the gluon's 4-momentum is not 
simply some multiple of that of the proton, the contribution from the photon-gluon fusion 
is a function of all four momentum components. In particular the kinematics of the virtual 
gluon production are important: where it is radiated from an on-shell parton, the gluon 
remembers the momentum fraction of the parton that emitted it. This, along with the 
discussion of the virtual gluon distribution in section |], indicates that we cannot detach 
the two halves of figure [l] but rather that the power corrections arise from the full diagram 
and are thus an 0(a 2 ) effect. 

Calculations of power corrections to structure functions and fragmentation functions due 



to gluon radiation by a quark are already published in the literature |8], ^, [H|. However 
the contributions from radiation by a gluon dominate, because of the relative behaviours 
of the quark and gluon parton densities. This larger correction was evaluated above for 
the structure functions. Above x = 0.05 the power corrections to the structure functions 
are small, but we predict that for x below 0.05 the corrections grow. So the correction to 
F2 at x = 0.01 and Q 2 = 4 GeV 2 might be around -2% (assuming a reasonable value of 
D\ m 0.1 GeV 2 ). The corrections to Ft and Fl are expected to be slightly larger. In any 
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case, we do not expect to see the large (~ 50%) corrections predicted by [§]. 

It was seen that the use of the pinch technique lead to a failure to recover the usual 
g — > g splitting function. This is not a problem — the remaining pieces come from the 
remaining diagrams. All the pinch technique has done is move certain terms from one 
diagram to another. But the point of using this particular separation of terms is that one 
can define a gauge-invariant QCD effective charge, so the virtuality of the gluon enters as 
the natural scale for the coupling; it is not clear that this is the case for the other terms. 

Finally, there is clearly some similarity between the situation here and that of a decaying 
outgoing virtual gluon in e + e~ annihilation (figure |5|), where we also have two renormalon 
chains with equal 4-momenta k, for some k. Again we have an integral over the virtuality 
k 2 , in which the coupling appears as \a s (— k 2 )\ 2 . But in that case we use equation ( |3.2| ) 
to convert a 2 into the spectral function p s (k 2 ), and the two chains became one, with a 
cut bubble insertion. There are two reasons why we cannot do the same here. Firstly, 
in the e + e~ case, the virtual gluon is timelike, k 2 is positive and thus we have p s (k 2 ) 
integrated over positive k 2 . But here the gluon is spacelike, so k 2 is negative, and while it 
is quite possible to convert the two factors of a s into a single p s {k 2 ), it is integrated over 
negative values of its argument. It is therefore not naturally manipulated into standard 
single-chain form. Secondly, although the algebra of the Feynman diagrams is identical 
(i.e. related by crossing), in the e + e~ case we can simply integrate over the cut bubble but 
in DIS we must include parton density functions for the incoming quark or gluon. This 
makes it quite impossible to integrate out the crossed cut bubble that is the mechanism for 
production of the virtual gluon. So, while there are interesting and useful parallels between 
power corrections to singlet DIS and e + e~ annihilation with outgoing gluon splitting, there 
are also significant differences and the relationship between them is not as simple as may 
naively have been supposed. 
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